ANALYTICAL RESULTS FOR 2-D NON-RECTILINEAR 
WAVEGUIDES BASED ON THE GREEN'S FUNCTION 
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Abstract. We consider the problem of wave propagation for a 2-D rectilinear optical waveguide 
which presents some perturbation. We construct a mathematical framework to study such a problem 
and prove the existence of a solution for the case of small imperfections. Our results are based on 
the knowledge of a Green's function for the rectilinear case. 
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1. Introduction. An optical waveguide is a dielectric structure which guides 
and confines an optical signal along a desired path. Probably, the best known example 
is the optical fiber, where the light signal is confined in a cylindrical structure. Optical 
waveguides are largely used in long distance communications, integrated optics and 
many other applications. 

In a rectilinear optical waveguide, the central region (the core) is surrounded by 
a layer with a lower index of refraction called cladding. A protective jacket covers the 
cladding. The difference between the indices of refraction of core and cladding makes 
possible to guide an optical signal and to confine its energy in proximity of the core. 

In recent years, the growing interest in optical integrated circuits stimulated the 
study of waveguides with different geometries. In fact, electromagnetic wave propaga- 
tion along perturbed waveguides is still continuing to be widely investigated because 
of its importance in the design of optical devices, such as couplers, tapers, gratings, 
bondings imperfections of structures and so on. 

In this paper we propose an analytical approach to the study of non-rectilinear 
waveguides. In particular, we will assume that the waveguide is a small perturbation 
of a rectilinear one and, in such a case, we prove a theorem which guarantees the 
existence of a solution. 

There are two relevant ways of modeling wave propagation in optical waveguides. 
In closed waveguides one considers a tubular neighbourhood of the core and imposes 
Dirichlet, Neumann or Robin conditions on its boundary (see |01j and references 
therein). The use of these boundary conditions is efficient but somewhat artificial, 
since it creates spurious waves reflected by the interface jacket-cladding. In this paper 
we will study open waveguides, i.e. we will assume that the cladding (or the jacket) 
extends to infinity. This choice provides a more accurate model to study the energy 
radiated outside the core (see |SL| and |Maj ). 

Thinking of an optical signal as a superposition of waves of different frequency 
(the modes), it is observed that in a rectilinear waveguide most of the energy provided 
by the source propagates as a finite number of such waves (the guided modes). The 
guided modes are mostly confined in the core; they decay exponentially transversally 
to the waveguide's axis and propagate along that axis without any significant loss 
of energy. The rest of the energy (the radiating energy) is made of radiation and 
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evanescent modes, according to their different behaviour along the waveguide's axis 
(see ^3] for further details). The electromagnetic field can be represented as a discrete 
sum of guided modes and a continuous sum of radiation and evanescent modes. 

As already mentioned, in this paper we shall present an analytical approach to 
the study of time harmonic wave propagation in perturbed 2-D optical waveguides. 
As a model equation, we will use the following Helmholtz equation (or reduced wave 
equation): 

Au(x, z) + k 2 n(x, z) 2 u(x, z) — fix, z), (1-1) 

with (x,z) £ K 2 , where n{x,z) is the index of refraction of the waveguide, k is the 
wavenumber and / is a function representing a source. The axis of the waveguide is 
assumed to be the z axis, while x denotes the transversal coordinate. 

Our work is strictly connected to the results in [MS a , where the authors derived 
a resolution formula for ( |1.1[ ), obtained as a superposition of guided, radiation and 
evanescent modes, in the case in which the function n is of the form 



n oi x ) = s , , (I- 2 ) 




where TLqq IS el bounded function decreasing along the positive direction and 2h is the 
width of the core. Such a choice of n corresponds to an index of refraction depending 
only on the transversal coordinate and, thus, ( |1.1[ ) describes the electromagnetic wave 
propagation in a rectilinear open waveguide. By using the approach proposed in |ACj . 



the results in |MS] have been generalized in |Cil] to the case in which the index 
of refraction is not necessarily decreasing along the positive direction. The use of 
a rigorous transform theory guarantees that the superposition of guided, radiation 
and evanescent modes is complete. Such results are recalled in S|3j The problem 
of studying the uniqueness of the obtained solution and its outgoing nature will be 
addressed elsewhere. 

In this paper we shall study small perturbations of rectilinear waveguides and 
present a mathematical framework which allows us to study the problem of wave 
propagation in perturbed waveguides. In particular, we shall assume that it is possible 
to find a diffeomorphism of R 2 such that the non-rectilinear waveguide is mapped in 
a rectilinear one. Thanks to our knowledge of a Green's function for the rectilinear 
case, we are able to prove the existence of a solution for small perturbations of 2-D 
rectilinear waveguides by using the contraction mapping theorem. 

In order to use such theorem, we shall prove that the inverse of the operator 



obtained by linearizing the problem is continuous (see Theorem 5.6 1. Such a problem 
has been solved by using weighted Sobolev spaces, which are commonly used when 
dealing with Helmholtz equation (see, for instance, |Le] ) . 

In a forthcoming work, the results obtained in this paper will be used to show 
several numerical results interesting for the applications. 

In f|2]we describe our mathematical framework for studying non-rectilinear waveg- 
uides. Since our results are based on the knowledge of a Green's function for rectilinear 
waveguides, in ^J3]we recall the main results obtained in [MS] , 

Section [3] will be devoted to some technical lemmas needed in ^5] The existence 
of a solution for the problem of perturbed waveguides will be proven in Theorem |5.6| 
Crucial to our construction are the estimates contained in f|5j in particular the ones 
in Theorem 15.41 
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AppendixfAjcontains results on the global regu larit y for solutions of the Hclmholtz 
; " r " N > 2, that we need in Theorem 



equation m 



5.6 



2. Non-rectilinear waveguides: framework description. When a rectilin- 
ear waveguide has some imperfection or the waveguide slightly bends from the recti- 
linear position, we cannot assume that its index of refraction n depends only on the 
transversal coordinate x. From the mathematical point of view, in this case, we shall 
study the Hclmholtz equation 



Au + k 2 n £ (x, z) 2 u = /, 



(2.1) 



where n s (x, z) is a perturbation of the function no(x) defined in ( 1.2 1, representing a 
"perfect" rectilinear configuration. 

We denote by L and L £ the Hclmholtz operators corresponding to n (x) and 
n £ (x,z) respectively: 

L = A + k 2 n (x) 2 , L £ = A + k 2 n £ (x,z) 2 . (2.2) 

In |MSj . the authors found a resolution formula for 

Lqu = f, 

i.e. they were able to write explicitly (in terms of a Green's function) the operator 
Lq 1 and then a solution of (1.1 1. Now, we want to use Lq 1 to write higher order 
approximations of solutions of (2.1|, i.e. of 



L £ u = f. 



(2.3) 



The existence of a solution of (2.3 1 will be proven in Theorem 5.6 by using a 
standard fixed point argument: since (2.3 1 is equivalent to 



Lqu = f + (L Q - L s )u, 



then we have 



u = L 1 / + eL Q 1 



U - L f 



Our goal is to find suitable function spaces on which Lq 1 and L °~ L ' are contin- 
uous; then, by choosing e sufficiently small, the existence of a solution will follow by 
the contraction mapping theorem. 

It is clear that this procedure can be extended to more general elliptic operators; 
in f|5]we will provide the details. 

3. A Green's function for rectilinear waveguides. In this section we recall 
the expression of the Green's formula obtained by Magnanini and Santosa in [MS 
and generalized in [Cil to a non-symmetric index of refraction. 

We look for solutions of the homogeneous equation associated to ( 1.1 ) in the form 

u(x,z) = v(x,X)e ikl3z ; 
v(x, A) satisfies the associated eigenvalue problem for v. 

v" + [A — q(x)]v = 0, inE, (3.1) 
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with 

n,=maxn, A = k 2 (n 2 — /3 2 ), q(x) = k 2 [n\ - n(x) 2 }. (3.2) 

IE 

The solutions of ([33J can be written in the following form 

{(j)j(h, A) cos Q(x — h) + ^ ] g' A ' > sin Q(x — h), if x > h, 
<j>i(x,\), ii\x\<h, (3.3) 

</)j(—h, A) cosQ(x + h) + ^^ A) sinQ(x + h), if x < -h, 



for j = s, a, with Q = >/X — d 2 , e? 2 = fc 2 (n 2 — n 2 z ) and where the <^-'s are solutions of 

(3.4) 



(3.1) in the interval (—h,h) and satisfy the following conditions: 

S (O,A) = 1, &(0,A) = 0, 



Q (O,A) = O, <^' a (0,A) = VA. 

The indices j = s,a correspond to symmetric and antisymmetric solutions, respec- 
tively. 



Remark 3.1. (Classification of solutions). The eigenvalue problem (3.1l leads to 
three different types of solutions of ( |1.1[ ) of the form up(x, z) = v(x, A)e 4fc/3z . 

• Guided modes: < A < or. It exists a finite number of eigenvalues X 3 m , 
m = 1 , . . . , Mj , satisfying the equations 

Vd 2 -Xcf> 3 (h,X) + ^(h,X) = 0, i6{s,a}, 



and corresponding eigenfunctions Vj(x,X J m ) which satisfy (3.1l. In this case, 
Vj(x, X J m ) decays exponentially for |a;| > h: 

'hiKX^e-^-^-V, x>h, 
Vj{x, X J m ) = <( (pj(x, X 3 m ), \x\ < h, 

4> j {~h,Xl n )e^ d2 - x3 "^ x+h \ x < -h. 

In the z direction, up is bounded and oscillatory, because /3 is real. 

• Radiation modes: d 2 < X < k 2 n 2 . In this case, up is bounded and oscillatory 
both in the x and z directions. 

• Evanescent modes: X > k 2 n 2 . The functions Vj are bounded and oscillatory. 
In this case (3 becomes imaginary and hence up decays exponentially in one 
direction along the z-axis and increases exponentially in the other one. 

By using the theory of Titchmarsh on eigenfunction expansions, we can write a 



Green's function for (1.1 ) as superposition of guided, radiation and evanescent modes: 

+ f°J\z-CWk*nl-\ 

G(x,z;t,Q= £ / Vj (x, X) vj (g, A) dpj (A) , (3.5) 



j<E{s,a} o 

with 

M 



" J 1 f v^A~ d 2 
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for all i] e C£°(R), where 

-j-oo 



Vj(x,\ J m ) 2 dx 



Vd 2 - XL 



Vd 2 - XL f ct> j {x,XL) 2 dx + ^{hiXin)' 



and where Vj(x, X) are defined by (3.3 1 (see |Cil] for further details). 
We notice that (3.51 can be split up into three summands 



G = G 9 + G r + G e , 



where 



G3(x,z;Z,0= E E 



M i e i\z-C\\/k 2 nl-\L 



je{sM}m=i 2i\J k 2 n 2 — XL 



j 

in • 



(3.6a) 



G r {x,z;i,0 = 



* J\z-<\y/k 2 nl-X 



2vr J 2iJk 2 n 2 - X 



■Vj(x, X)vj(£, X)aj(X)dX, (3.6b) 



+00 

G e (x,z;£,0 = -^ E / : 



,-\z-(\y/\-k2nl 



J 2^/X-k 2 n 2 



with 



(A-^)^(/i,A) 2 + ^(M) 5 



/(a;, A)u,-(£, A)aj(A)dA, (3.6c) 



(3.7) 



G 9 represents the guided part of the Green's function, which describes the guided 
modes, i.e. the modes propagating mainly inside the core; G r and G e are the parts 
of the Green's function corresponding to the radiation and evanescent modes, respec- 
tively. The radiation and evanescent components altogether form the radiating part 
G rad of G: 



G rad = G r + G e = 



E 



1 



2n -Jr 1 x J 2iJk 2 nl - A 



v J (x,X)v j (^,X)a j (X)dX. (3.8) 



4. Asymptotic Lemmas. This section contains some lemmas which will be 
useful in the rest of the paper. 

Lemma 4.1. Let q G Lj oc (M) and Xo = min(Af,Af), where Af and A" are defined 
in Remark 3.1 Let <f>j, j € {s,a}, be defined by (3.3 I. Then, the following estimates 
hold for x g \—h, h] and A > Xq: 

|&(*,A)| < 



i>'Jx,X)\ < <S>«Vx, 



(4.1) 
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whe 



<f>» := exp 



\q(t)\dt 



(4.2) 



Proof. We consider the function 



and notice that 



i>j(x, A) = (jy'Jx, A) 2 + \(j)j{x, \y 



ip'Ax, A) = 2q(x)(j) j (x, X)(j)'Ax, A), 



as it follows from (3.1). By using Young's inequality, we get that ijjj satisfies 

^(x,X)<^^(x,X), 
>,(0,A) = A. 

Therefore, by integrating the above inequality, we obtain that 

h 

A) < A exp { ^= I \q{t)\dt } < A$ 2 



'X 



which implies (4.1|. □ 

In the next two lemmas we study the asymptotic behaviour of the function Oj (A) 
as A — ► +oo and A — * d 2 , respectively. 



Lemma 4.2. Let <7j(A), j € {s,a}, be the quantities defined in (3.7l. The folloi 



ing asymptotic expansions hold as X — > oo: 



1 



, a a (X) = - f +(D 
V A 



(4.3) 



Proof. By multiplying 

<t>"(x, A) + [A - q{x)](j)j{x, X) = 0, x e [-ft, ft], 
by <^-(x, A) and integrating in a: over (0,/i), we find 

^(/ 1 ,A) 2 -^.(O,A) 2 + (A-d 2 )[0 3 (/ l ,A) 2 -0 J (O,A) 2 ] -2 / [g(a:)-d%(a:, A)^(x, A)dx 



Thus, by using (3.4 1, we obtain the following inequalities: 



U' s (h, A) 2 + (A - d 2 )cf> s (h, A) 2 - (A — d ) < 2fc 2 (< + n 2 ,) / |^(x, A)#(a:, A)|das 



U;(^,A) 2 + (A-d 2 )0 Q (/ l ,A) 2 - A| <2fc 2 (7i 2 + n 2 ; ) / \4> a (x, X)(/>' a (x, X)\dx. 
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The asymptotic formulas ( |4.3[ ) follow from the two inequalities above, (3.7 1 and 
the bounds (4.1) for <fij(x,A) and (j>'Jx,A). □ 



Lemma 4.3. Let crj(A), j £ {s,a}, be the quantities defined in ( 3.7 1 . The follow- 
ing formulas hold for A — > d 2 : 



a,- (A) = < 



' V^d 2 



+ O (A - d 2 



if<t>'Jh,d 2 )^0, 



(4.4) 



— d 2 ^j , otht 



>j(M 2 ) 2 \/A~^2 

Proof. We recall that, if q £ Lj oc (R), for x £ [— /i, /i], (f) s (x, A) and a (x, A) are 
analytic in A and vA, respectively (see |CLj ). Thus, in a neighbourhood of A = d 2 , 
we write 

+oo +oo 

4>(h, A) = £ (A - rf 2 ) m a m , ^(fc, A) = £ (A - d 2 ) m b m , 



m=0 



4>{h,X) 2 = £(A-d 2 ) m a m , 0'(M) 2 = £(A-d 2 ) m A 



m=0 



where we omitted the dependence on j to avoid too heavy notations. 



We notice that o.q = a 2 ,, ol\ — 2ao&i and the same for bo and b\. From (3.7) we 



have 



aj(X)~ 1 = ^/\-d 2 ^(h,\) 2 + 



A, 



V>^d 2 J 



(4.5) 



m— 



If Po 7^ 0, since /3o — Oq, (4.4) follows. If /3q = we have that the leading term in 
(4.5) is ccq + P\. We notice that ao ^ 0, otherwise (j)(x,d 2 ) = for all x £ K. We 
know that ^ = 2b bi = 0, because 6 = 0- Then a + fi\ = a and (14.4b follows. □ 



5. Existence of a solution. Let /i : M 2 — > R be a positive function. We will 
denote by L 2 (n) the weighted space consisting of all the complex valued measurable 
functions u(x, z], {x,z) £ M 2 , such that 

yflu £ L 2 (R 2 ), 

equipped with the natural norm 

IMIi*0*) = / \u(x,z)\ 2 fi(x lZ )dxdz. 

R 2 

In a similar way we define the weighted Sobolev spaces H 1 ^) and H 2 {ff). The 
norms in H 1 ^) and H 2 {^jl) are given respectively by: 



IMIi^Gu) = / l"^ 2 )! f-( x i z)dxdz + / |Vw(a;,z)| (i(x,z)dxdz 
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and 

II w IIh 2 (/j) = J \u(x, z)\ 2 n(x, z)dxdz+ J \Vu(x, z)\ 2 /i(x, z)dxdz+ J |V 2 it(x, z)\ 2 [i(x, z)dxdz. 

R 2 R 2 R 2 

Here, Vw and V 2 u denote the gradient and Hessian matrix of u, respectively. 

In this section we shall prove an existence theorem for the solutions of (2.3 1. We 
will make use of results on global regularity of the solution of (1.1); such results will 
be proven in Appendix [A") 

The proofs in this section and in Appendix [A] hold true whenever the (positive) 
weight (i has the following properties: 



|V M | < C lt i, |VV| < Cm, in M 2 , 

where C\ and C 2 are positive constants. 

In this section, for the sake of simplicity, we will assume that /1 is given by 



(5.1) 



li{x,z) = ni{x)ii 2 {z)- 



(5.2) 



with fij e X°°(M) n L 1 (M), j = 1,2. Analogous results hold for every [i satisfying 
(5.1 ) and such that 

fji(x, z) < ix 1 {x)n 2 {z). (5.3) 

For instance, it is easy to verify that the more commonly used weight function 
n(x, z) = (1 + x 2 + z 2 )~ a , a > 1, satisfies (5.1 1 and (5.3 1. 

Before starting with the estimates on u, we prove a preliminary result on the 
boundness of the guided and radiated parts of the Green's function. 



Lemma 5.1. Let G 9 and G r be the functions defined in (3.6al and (3.6b), re- 
spectively. Then 



Mi 



\G 9 (x,z;^C)\<^l E E 



je{s,a}m=i 2\ k 2 n1 — A 



(5.4a) 



and 



\G r (x, z; £X)\ ^ max ' 



y r Jt — y t 2 

2' 4J^r ^ 37 2tt ^ 3 
je{s,a} je{s,a} 



(5.4b) 



Here, 



1 * „ 
K 1 1 „ 



a, (A) 



2s/k 2 nl - A 



dX 



j e {s,a}, 



where, as in Lemma 4-L is given by (4.2). 

Proof. Since G 9 is a finite sum, from Remark |3.1| and Lemma |4.1[ it is easy to 
deduce (5.4a I. 
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In the study of G r we have to distinguish three different cases, according to 
whether (a;, z) and (£, Q) belong to the core or not. Furthermore, we observe that 



Ti < +00 as follows form Lemma 4.3 



Case 1 : x,£ € [—h,h]. From Lemma 4.1 we have that u, (x, A) are bounded by 
<f>*. From ( 3.6b I we have 

\G r \<l E 

je{s,a} 



Case 2 : |x|, |£| > h. We can use the explicit formula for Vj (see (3.3)) and obtain 
by Holder inequality 

\vj( x> \)\ < ^/mm^ + Q-V/m) 2 = [Qo-i(A)]-*. 

Therefore we have 



\G r (x,z;{,0\<± £ J 

j€{s,a} p 



dX 



\J 2VX^Jk 2 nl-X 2' 

a} d 2 V 



and hence (5.4b I follows. 



Case 3 : |x| < h and |£| > fa. We estimate \vj(x, X)\ by <I>* and Uj(£,A) by 



[Qcrj(A)] 2 ) and write: 



\G r (x,z;t,()\<-<S>* 

zir * — • 



j'efs 



a} d 2 V 



Qj(A) 

Q 



dX 

VX^&y/W^~X 



Again (5.4b I follows. □ 



In the next lemma we prove estimates that will be useful in Theorem |5.4| 

Lemma 5.2. Let m € L 1 (R) and ^2 € ^(R) flL 2 (K). For jj € {s,a} and 
A, 77 > k 2 n 2 , set 



+00 



pj,i(X,rj) = / Vj(x,\)vi(x,rj)fii(x)dx, 



(5.5) 



and 



+00 



q(X,7}) = / (eA,,, -k ^2){z)^ 2 {z)dz 1 

— OO 

where e X}V {z) = e -M(\A 3 ^+\A r ^) . 



(5.6) 
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Then p jt i(X, r?) = for j ± I, 



x/X^d 2 



!, (5.7) 



and 



\q(X,n)\ < min \\n2Wl, 



</A - fc 2 n 2 ^ _ fc2 n 2 



(5.8) 



Proof. Since v s and u a are, respectively, even and odd functions of x, then 



By Holder inequality, pjj(X,r]) 2 < Pj,j(X, \)pjj(r],r]); also, by the formula (3.3) 
we obtain 

h +00 

\p j , j (X,X)\<2 [ \4> j (x,X)\ 2 \» 1 {x)\dx + 2 [ |« i (as,A)| 2 |/ii(as)|da; 



-Hoc. 



< 2$ 2 / \ ftl (x)\dx + 



VX^&a^X) 



\H\(x)\dx, 



and hence (5.7) 



Now we have to estimate q(X,rf). Firstly we observe that 



+00 



+00 



\q(X,v)\ < / ^{z)dz / M2(CX=llM2||? 



which proves part of (5.8 1. Furthermore, from Young's inequality (see Theorem 4.2 
in |LLp and the arithmetic-geometric mean inequality, we have 



|?(A,77)| < ||e A ^||i||/i2|| 



2|M|1 



< 



v/A — fc 2 n 2 + \Jrj — k 2 ri 
yA — fc 2 n 2 \]r\ — k 2 n 2 ' 



which completes the proof. □ 



Theorem 5.3. Let G be the Green's function (3.51. TTien 

\\G\\L 2 (fiXfi) < +°°- 



(5.9) 



Proof. We write G = G 9 + G r + G e , as in ( |3.6a| )-( |3.6c) , and use Minkowski 
inequality: 



\\G\\l 2 (^x^) < I|G 9 || L 2 (aiXai) + ||G r ||_ £ ,2( /JX(tt) + ||G e || L 2 (MXAl) . 



(5.10) 
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From Lemma 5.1 and (5.1) it follows that 



\\G 9 \\l 2 (pxh), \\G r \\L2(px fJ ,) < +oo- 



(5.11) 



It remains to prove that 1 1 C e 1 1 ^2 x ^j,) < +oo. From (5.2) we have 
l|G e || 2 L2M = / \G%x,z;Z,0\ 2 »i(x)M(z)»i(0»2(0dxdzdtdt 



G e (x, z; f, ()G e (x, z; f , ()/j,i(x)[x 2 (z)lJ,i(£,)lJ,2(0 dxdzd ( d C, 



hence, thanks to Lemma 5.2 the definition (3.6c I of G e and Fubini's theorem, we 
obtain: 



— y 



+00 +00 



q(X, rj) dXdrj 



\/A — fc 2 n 2 W rj — k 2 n 2 



The conclusion follows from Lemmas 5.2 4.2 and 4.3 □ 



Corollary 5.4. Let u be the solution of (1.1 1 given by 



u(x,z)= G(x,z;(;,(;)f(£,()d£d(, (x,z) € 



(5.12) 



with G as in (3.5 1 and let f £ L 2 (/j, 1 ). Then 

IMIff^) < C'H/IU^-i), 

where 



C 2 = - + 2C 2 



+ 4C 2 + 8C| + (1 + 4C 2 )fc 2 n 2 + 2fc 4 n 4 



|G| 12 



(5.13) 



(5.14) 



Proof. From (5.121 and by using Holder inequality, it follows that 



\W\\h 



G(x,z;£,()f(£X) d £ d ( fJ-(x,z)dxdz 



< / / |G(x,z;^C)| 2 M(e,C)M^^)^C^ 

R 2 R 2 

= II/IIl 2 ( AI - i )II^'IIl 2 ( A1 x AI )' 



and then we obtain ( |5.13 1 and (5.14) from Lemmas A.l and A. 3 □ 



Remark 5.5. In the next theorem we shall prove the existence of a solution of 
L £ u = f. It will be useful to assume in general that L £ is of the form 



2 



(5.15) 



i=i 
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This choice of L £ is motivated by our project to treat non-rectilinear waveguides. 
Our idea is that of transforming a non-rectilinear waveguide into a rectilinear one by 
a change of variables r : M 2 — ► M 2 . 

For this reason, we suppose that T is a C 2 invertible function: 

T(s,t) = (x(s,t),z(s,t)). 

By setting w(s, t) — u(x, z), a solution u of ( |1.1[ ) is converted into a solution w of 

|Vs| 2 w ss + | Vt\ 2 w tt + 2Vs • Vt w st + As ■ w s + At ■ w t + c(s, t) 2 w = F(s, t), (5.16) 

where c(s, t) = kn(x(s, t), z{s, t)) and F(s, t) = f(x(s, t), z(s, t)). 

If our waveguide is a slight perturbation of a rectilinear one, we may choose T as 
a perturbation of the identity map, 

F(s,t) = (s + e<p(s,t),t + eip(s,t)), 

and obtain L e w = F from ( |5.16 I, where L e is given by (5.15), with 

afj = Sij + ea^j, i,j = 1,2; b, = e&f, i = 1,2; c e = /c 2 n (a;) 2 + ec e ; (5.17) 

we also may assume that 



E 



i=l 



|c e | < in 



for some constant K independent of e. 



(5.18) 



Theorem 5.6. Let L £ be as in Remark 5.5 and let f G L 2 (/i _1 ). T/ien £/iere 
exists a positive number Eq such that, for every e G (0,£q), equation L £ u = f admits 
a (weak) solution u e G if 2 (/i). 

Proof. We write 

L e = L + eL e ; (5.19) 

clearly, the coefficients of L e are of-, 6f and c e defined in (5.171. We can write ( 5.19| ) 
as 

M + £Lfi 1 L e u = L„ 1 f; 

L^ x f is nothing else than the solution of (1.1) defined in (5.121. 

We shall prove that UQ X L e maps H 2 (n) continuously into itself. In fact, for 
u G H 2 (fi), we easily have: 



\L F u\ 



E (^) 2 E 



Lt,j=l 



»J=1 



2 2 

i 2 +E^) 2 En 2 + ( ££ ) 2 m 2 

i=l i=l 



/i 1 dxdz 



< K 2 \\u\\ 2 H2M . 



Moreover, Corollary 5.4 implies that 

ll-^o /lli? 2 (ju) - ^ll/IU 2 ^- 1 )' 
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and hence 

ll-^o 1 Leu\\h*(p) < CK\\u\\ H 2 (m) . 

Therefore, we choose e = (CK)^ 1 so that, for e S (0,e ), the operator sLq 1 L e is a 
contraction and hence our conclusion follows from Picard's fixed point theorem. □ 



6. Numerical results. In this section we show how to apply our results to 
compute the first order approximation of the solution of the perturbed problem. The 
example presented here has only an illustrative scope; a more extensive and rigorous 
description of the computational issues can be found in |Cil| and |Ci2j . where we 
apply our results to real-life optical devices. 




(a) The perturbed waveguide. The dashed lines (b) Real part of w^°\ Here, ui' ' is a pure 

show the effect of T on the plane. In particular, guided mode supported by the waveguide in the 

they show how a rectangular grid in the (s, t)- rectilinear configuration, 
plane is mapped in the (x, z)-plane. 



Figure 6.1. The perturbed waveguide and the real part of w'°> . 



In this section we study a perturbed slab waveguide as the one shown in Fig 6.1(a) 



In the case of a rectilinear slab waveguide it is possible to write the Green's formula 
explicitly and numerically evaluate it (see |MS| ). 

Having in mind the approach proposed in Remark |5.5[ we change the variables 
by using a C 2 — function r : R 2 — > M 2 of the following form: 

T(s,t) = (s,t + eS(s)T(t)), 



where S,T S C 2 (R); a good choice of S and T is represented in Fig 6.2 In Fig 6.1(a) 



we also show how T transforms the plane, by plotting in the (x, z) — plane the image 
of a rectangular grid in the (s, t)— plane. 

By expanding L E and w by their Neumann series, we find that and w^ 1 ' (the 
zeroth and first order approximations of w, respectively) satisfy 

Aw^ + k 2 n{s) 2 w^ = F(s, t), (6.1a) 

and 

AiijW + k 2 n(s) 2 w^ 

= -2S'(s)T(t)wi°J - 2S( S )r(t)wf t ) - [S"(s)T(t) + 5( S )r"(t)] Wt (0) , (6.1b) 
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25'h S-aXIS 



(a) The function S. 



(b) The function T. 



Figure 6.2. Our choice of the functions S and T. Such a choice corresponds to a perturbed 
waveguide as in Fig \6. 1 (a)\ 





(a) Real part of u>W. 



(b) Modulus oftuf 1 '. 



Figure 6.3. The real part and modulus of «)W (the first order approximation of us) 



respectively. 

In our simulations, we assume that is a pure guided mode and calculate 



by using (6.1 1 and the Green's function (3.5 1. In other words, we are taking a special 



choice of / and see what happens to the propagation of a pure guided mode in the 
presence of an imperfection of the waveguide. 

In Figures [6~I(b)1 16~3| and |6~4| we set k = 5.0, h = 0.2, n co = 2, n d = 1. With 
such parameters, the waveguide supports two guided modes, corresponding to the 
following values of the parameter A: Xf = 23.7 and A° = 73.5. 

As already mentioned, we are assuming that is a pure guided mode. Here, 



-,(o) ; 



is forward propagating and corresponds to A^: 



the real par t of is s hown i n Fig 6.1(b) 
and 



Figures 6.3(a) 



6.3(b) show the real part and the absolute value of wr 



respectively. We do not write here the numerical details of our computation and refer 
to |Cil] and |Ci2j for a more detailed description. 

In Figures 6.4(a)| and 6.4(b)| we show the real part and the absolute value of 
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(a) Real part of io(°) + ewW. (b) Modulus of u>(°) + ewW. 



Figure 6.4. The real part and modulus of iu(°) + ew^ . TTie pictures clearly show the effect 
of a perturbation of the waveguide: due to the presence of an imperfection, the waveguide does not 
support the pure guided mode w( ' and the other supported guided mode and the radiating energy 
appear. 

w^ ' + ew^\ respectively. Here, we choose e = 1 to emphasize the effect of the 
perturbation on the wave propagation. As is clear from Theorem |5.6| our existence 
result holds for e £ [0,£o], where Eq = (CRT) -1 (which will be presumably less than 
1). The computation of £o an d the convergence of the Neumann series related to w 
have not been considered here; again, we refer to |Cil] and |Ci2] for a detailed study 
of such issues. 

7. Conclusions. In this paper, we studied the electromagnetic wave propagation 
for non-rectilinear waveguides, assuming that the waveguide is a small perturbation 
of a rectilinear one. Thanks to the knowledge of a Green's function for the rectilinear 
configuration, we provided a mathematical framework by which the existence of a 
solution for the scalar 2-D Hclmholtz equation in the perturbed case is proven. Our 
work is based on careful estimates in suitable weighted Sobolev spaces which allow us 
to use a standard fix-point argument. 

For the case of a slab waveguide (piecewise constant indices of refraction) , numer- 
ical examples were also presented. We showed that our approach provide a method 
for evaluating how imperfections of the waveguide affect the wave propagation of a 
pure guided mode. 

In a forthcoming paper, we will address the computational issues arising from the 
design of optical devices. 



Appendix A. Regularity results. In this section we study the global regularity 
of weak solutions of the Helmholtz equation. Since our results hold in R*, N > 2, it 
will be useful to denote a point in by x, i.e. x = {x\, . . . , Xn) £ M. N . 

The results in this section can be found in literature in a more general context 
for N > 3 (see |Ag| ) . Here, under stronger assumptions on n and /i, we provide an 
ad hoc treatment that holds for N > 2. 

We will suppose / £ L 2 (^ _1 ). Since fi is bounded, it is clear that / £ L 2 (/j.) too. 
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Lemma A.l. Let u 6 Hf oc (R N ) be a weak solution of 
Au + k 2 n(x) 2 u = f, xeR N , 
with n £ L°°(R ). Let /i satisfy the assumptions in (5.1|. Then 

j '\Vu\ 2 ndx<~ J \f\ 2 fidx+^2C 2 + k 2 nl + ^j J \u\ 2 fidx, 



(A.l) 



(A.2) 



where n* = ||n|| £ =o( R iv). 

Proof. Let 77 G Cf?(R N ) be such that 

77(0) = 1, < 77 < 1, |Vt7|<1, |V 2 7 ? |<1, 

and consider the function defined by 

H m {x) = fjb(x)ri (—) ■ 
Then /x m (x) increases with to and converges to n(x) as m — > +00; furthermore 



(A.3) 



(A.4) 



\V(* m (x)\ < |Vm(»)| + — A»(a?) < ( Ci + — J /a(x), 
to \ m J 

|V 2 M m(x)| < |V a A*(a;)| + -|V M (z)| + -A^fr) < (c 2 + — + \ ) 



TO TO 



(A.5) 



for every a; G 



Since u is a weak solution of (A.l ), we have that 



fvu-V^dx-k 2 Jn(x) 2 uHx=-J 

I" 1" R N 



f<t>dx, 



for every (j) 6 iJj^QR^). We choose (/> = and obtain by Theorem 6.16 in [LLj : 
J \Vu\ 2 /i m dx — — J uVu ■ V [i m dx + k 2 J n{x) 2 \u\ 2 pi m dx ~ J fu/i m dx. (A. 6) 



R N 



Integration by parts gives: 

Re / uVu • S//i m dx = 



~ J |u| 2 A/j m da;; 



hence, by considering the real part of (|A.6|), we obtain 
j \Vu\ 2 ii m dx < 2 (c 2 



2C X 1 



TO 777/ 



\u\ /idx + k / |u| fJ, m dx 



R N 



\f\ 2 H m dx 



\u\ 2 ii m dx 
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here we have used (A. 5) and Holder inequality. 



Young inequality and the fact that \x m < /i then yield: 



|V u\ 2 \x m dx < 



2 [Co 



2Ci 1 



m m 



2„2 



k n 



\u\ 2 fidx + ^ / \f\ 2 fidx. 



The conclusion then follows by the monotone convergence theorem. □ 

Lemma A. 2. The following identity holds for every u £ Hf oc (M. N ) and every 
<f>£C 2 (R N ): 



J \Au\ 2 (f>dx + J \\7u\ 2 A<j)dx = J \\7 2 u\ 2 cj)dx + Re J (\7 2 <f)X7u,X7u)dx. (A.7) 



Proof. It is obvious that, without loss of generality, we can assume that u £ 
C^K^); a standard approximation argument will then lead to the conclusion. 

For u £ C 3 (R N ), (A.7l follows by integrating over M. N the differential identity 



UiUj 



N N N N 

uuujj o u a^ij + X! "i^fe - Re 

*,J=1 ij'=l *,J=1 i,j=l 



and by divergence theorem. □ 



Lemma A. 3. Let u £ H 1 {^jl) be a weak solution of (A.l). Then 



J \V 2 u\ 2 ixdx < 2 J \f\ 2 fidx + 2k A nt J \u\ 2 ^dx + 4C 2 J \Vu\ 2 ^dx. (A.8) 

K N R N R N R N 



where C2 is the constant in (A. 5) 



Proof. From well-known interior regularity results on elliptic equations (see The- 
orem 8.8 in [GTj ) . we have that if u £ H\JM. N ) is a weak solution of (A.l I, then 



u £ Hf oc (R ). Then we can apply Lemma 



A.2 



to u by choosing <p = /i J( 



J \Au\ 2 /j, m dx+ J \\7u\ 2 A/i m dx — J \V 2 u\ 2 ^, m dx + Re J (V 2 /i m Vu, Vu)dx. 

R N K N R N E N 
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From (A.l I, (A. 5 1 and the above formula, we have 



J \\7 2 u\ 2 n m dx = J |/ - k 2 n(x) 2 u\ 2 ^ m dx + J \\7u\ 2 A^, m dx - 2Re J (S7 2 n m \7u,S7u)dx 
<lj \f\ 2 n m dx + 2k i nt J \u\ 2 [i m dx 



1" R N 



J \Vu\ 2 Lidx + 2 J |VVm||Vw| 2 efo 



<2 / \f\ 2 f i m dx + 2k 4 nt I \u\ 2 ti m dx 

E 

m m z 

Since fi m < fi, the proof is completed by taking the limit as m — > oo. □ 
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